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We study a three dimensional Abelian Higgs model containing singly- and doubly-charged scalar 
fields coupled to a compact Abelian gauge field in the London limit. The model attracts interest 
because of its relevance to high-r c . superconductors with charge 1 holon and charge 2 spinon-pair 
fields. It contains two types of vortices carrying magnetic flux and one type of instanton-like 
monopoles. Using thermodynamic and topological observables we present the phase diagram in 
the parameter space of the gauge and holon and spinon-pair couplings. The Fermi liquid, the spin 
gap, the superconductor and the strange metallic phases have been identified in a wide region of 
parameters. The model may serve as a toy system modelling non-perturbative properties of the 
Yang-Mills theory. 
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1. Introduction 

The physics of high- T c superconductivity [|I|] is not understood yet. At normal temperatures, all 
known high-r r superconductors are ceramic crystals characterized by a poor conductivity. At low 
temperatures the clean ceramic materials are rather insulators than conductors (Mott insulators). 
However, as one adds impurities to the clean material ("doping"), a good insulator becomes in 
certain cases a superconductor at low enough temperatures. 

Despite superconducting specimen are three-dimensional structures the physics of high-r c su- 
perconductivity is believed to be essentially two-dimensional [JTJ] _ In fact, all known high-r f super- 
conductors consist of copper oxygen (CUO2) lattice planes. In the undoped state the crystal does 
not contain enough free carriers of electric charge. The role of impurities is to provide the carriers 
- electrons (like in Nd2-. v Ce v Cu04 case) or holes (like in La2- x Sr t Cu04 material) - into the CUO2 
planes, eventually making the CUO2 planes superconducting 2D systems. Therefore, the physics 
of high-T c superconductivity must be understood as an essentially 2D phenomenon. 

Physically, there are two basic parameters: the temperature T and the concentration of im- 
purities (doping) x. The phase diagram of a real high-r c . superconductor is the temperature- 
concentration plane and typically contains the following phases: a metallic, a pseudogap (found 
in hole doped materials), an anti-ferromagnetic and a superconducting phase. 

One popular approach to the physics of superconducting planes is based on a simplified model 
describing the dynamics of holes and spins. Basically we have a lattice with hopping holes and lo- 
calized spins acting as dynamical variables. The lattice spacing and possible anisotropy is dictated 
by the ceramic host material itself. The effective description is provided by the t — J Hamilto- 
nian [§: 

H t j = -t Y, c L( l - n h-o){l-rij-o)cj0+J'£{S i Sj--n i nj), (1.1) 

<ij>,a <ij> 

where the first term describes the hopping of holes or electrons without changing spin while the 
second term describes the anti-ferromagnetic Heisenberg coupling between the spins located at the 
copper-sites. Here S,- = (l/2)£ ffff /cJ ff cr ffCT 'C; ff ' is the spin operator, c ;(T and «; jf7 are the electron 
creation and occupation number operators for given spin a, respectively, and ?i, = Lcj c L c 'ff- 



Even the simplified t — J model fll.ip is difficult to solve. One successful approach is based 
on the slave boson formulation which proposes to split the spin and charge variables of the 
electrons. Let the electron creation operators be written as c] a = fj a bj, where fi a is a spin-particle 
("spinon") operator and b[ is a charge-particle ("holon") operator. In order to forbid double oc- 
cupancy of sites one imposes the constraint ftft\ +fj\fii + b]b{ = 1 on the physical states of the 
system. 

The spinon is a fermionic chargeless particle which carries information about the spin of the 
electron while the holon is a bosonic particle which is responsible for the electron charge. The 
essential feature of the spin-charge separation approach is that it introduces an additional (internal) 
compact U (1) degree of freedom which can be formulated as a gauge freedom on the operator 
level: 

c ia -> e' r, c ia : fia ^ fia, bi^e' y, b ia [usual gauge freedom] , (1.2) 
eta ->■ c ia : fia -> e m f ia , bi^e ia, b ia [internal gauge freedom] . (1.3) 
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2. The compact (7(1) gauge model coupled to matter 

The emerging effective theory of superconductivity can further be simplified and reformulated 
as a lattice gauge model [0]. The basic idea is to neglect - or absorb into the couplings of the 
effective model - the external gauge degree of freedom (|1.2[) and to concentrate attention to the 



internal compact gauge degree of freedom ( |1.3| ). In fact, the coupling of the usual electromagnetic 
field to the charge-carrying holon degrees of freedom (as well as the inter-holon interaction due to 
the electromagnetic interaction) is relatively weak compared to the strong correlations among the 
electrons in the lattice environment. 

Thus, we go over from the t — J model ( |1.1| ) to a compact Abelian gauge model with internal 



symmetry ( |1.3| ), which couples holons and spinons. As in usual BCS superconductivity, at certain 
parameters of the t — J model the spinons couple and form bosonic quasiparticles. In a mean field 
theory one can define the fields 



Xij — ^(/iff/ycr) "~ * Xij ' e 



-i(ai-aj) 



{fnfii- fnfii) 



Jidj+Clj) 



(2.1) 



which behave under the internal gauge transformations (1.3) as a neutral (vector-like) particle and 
a doubly-charged matter field, respectively. The phase of the # -particle is associated with the 
compact U (1) gauge field, 6{j = axgXij ~ * % — a i + a j- The radial part of the # -particle defines 
the so-called "resonating valence bond" (RVB) coupling % = (\Xij\), and the doubly-charged matter 
field A is called spinon-pair field (an analog of the doubly-charged Cooper pairs). 

Usually the RVB coupling is treated in the mean-field approximation and therefore is assumed 
to be fixed, %(x) = % = const. Thus, the dynamical content of the effective model is given by a 
singly-charged boson field (holon) b, a doubly-charged boson field (spinon-pair) A and the compact 
(7(1) gauge field d. 

At high temperature the RVB coupling is vanishing, % = 0, and the system is in the Mott 
insulator (or, poor metallic) phase, see Fig. |l[ As the temperature decreases, the RVB coupling is 
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Figure 1: The schematic phase diagram of the t/(l) model for a high-r c - superconductor 



getting non-zero, % / 0, enabling the formation of the spinon-pair condensate A = |(Ay)| and/or 
of the holon condensate b = (bi). Depending on the presence of the condensates, four phases 
classified in Ref. [^] may emerge. In Fig. [l] they are sketched and denoted as the Fermi liquid, the 
spin gap, the superconductor, and the strange metallic phases. 
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3. The compact Abelian two-Higgs model on the lattice 

Following the proposal in Ref. ^ the model described in Section || can be studied as compact 
Abelian two-Higgs model (cA2HM) in three dimensions with a U{\) gauge link field 0/, a single- 
charged holon field <&i, and a double-charged spinon-pair field <I>2. The qualitative characteristics 
of the model become already visible in the London limit, in which the radial parts of both Higgs 
fields <!>; = l^ le"^ are frozen, |^>,-| = const, i = 1,2. The remaining dynamical Higgs variables are 
the phases of the matter fields (pi^. The action of the cA2HM model in the London limit is: 

S[d, (pi ,(pi} = -P £cos e P - K X £cos(d(p 1 + 0), - K- 2 £cos(d< ? > 2 + 20)/ , (3.1) 
p I I 

where dp is the standard lattice plaquette field, d denotes the ordinary lattice derivative, (d<p) v ^ = 
(Px+fi — <Px- The parameter /3 is the inverse gauge coupling, fq oc t ■ x and ^ °< 7 are the hop- 
ping parameters of the holon and the spinon-pair field, respectively. They are proportional to the 



parameters of t — J model. The model (3.1) obeys the U{\) gauge invariance: 



0^0 + da, <pi^<pi-a, <p2^<P2-2a, (3.2) 



which is a direct counterpart of the internal gauge symmetry (|1.3|). 

In the limit K"2 — ► one gets the compact U{\) gauge model with a charge- 1 Higgs field only, 
the holon. The phase diagram of the 2=1 compact Abelian Higgs model contains [|], ||] two 
regions: the confinement "phase" (characterized by a low holon condensate b) at small hopping 
parameter Ki) and the Higgs "phase" at large K\ (where the condensate b becomes large). Both 
regions are analytically connected for strong enough Higgs selfcoupling, although a Kertesz line [||] 
physically separates them also there. A qualitatively similar picture (however with a true second 
order phase transition) is realized [^] in the limit K\ — > where one gets the Q = 2 compact Abelian 
Higgs model in which the spinon-pair field plays the role of the sole Higgs field. Thus, we expect 
that the phase diagram of the model (3T) should contain all four phases depicted in Fig. [T] in the 
X > region. 

Below we report our Monte Carlo investigation of the cA2HM which is still under way. To 
scan the phase diagram we have simulated on a 16 3 lattice, choosing three values of the gauge 
couplings, j6 = 1.0, 1.5 and 2.0, for a huge grid of Ki^ hopping parameter pairs covering the 
range 1 < K"i 2 ^ 2.5. We are particularly interested in the strongly coupled case, j8 ~ 1. 

The compactness of the gauge field guarantees the presence of instanton-like magnetic mono- 
poles. An elementary monopole is a source of F mon = 2% units of magnetic flux associated with the 
internal gauge field 0. Due to the matter fields two types of topologically stable vortices exist. The 
magnetic flux quanta of the vortices corresponding to the holon and the spinon-pair "Higgs" fields 
are f , ™ rt = 2% and F 2 V01t = n, respectively. Since the magnetic flux is conserved, one monopole is 
simultaneously a source of one holon vortex and two spinon-pair vortices. 

The densities of the three topological defects are logarithmically plotted in Fig. ^| over the 
(k"i , K2) plane for a strong gauge coupling j8 = 1. With increasing hopping parameters the monopole 
density gets suppressed. The density of the holon (spinon-pair) vortex becomes suppressed with 
increasing holon (spinon-pair) hopping parameter K\ (K2), beyond a line in the Kj-K^ plane. 



'The actually relevant region of the ffi-Kj plane depends on jS. 
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monopole holon vortex spinon-pair vortex 

Figure 2: The densities of the topological defects over the (ki, K2) plane at /3 = 1 . 



The connectivity of the vortex clusters gives a clear view of the phase diagram of the model [g, 
|7[ |6p. A Higgs condensate suppresses the proliferation of the corresponding vortices: they are pre- 
vented to percolate over infinitely long distances. We show in Fig. |3| the percolation probabilities 
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Figure 3: The percolation probability of the infrared vortices and the phase diagram of the system at /3 = 1 . 

dropping to zero (extracted from the cluster correlation functions) and the corresponding phase 
diagram (classified according to Ref. [gj, see Fig. [I]). To identify the phases, we measured the 



average plaquette and both link contributions of the action and their susceptibilities ( p.l| ), respec- 
tively. Some examples are shown in Fig RL From a preliminary finite size analysis from 12 3 to 32 3 
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Figure 4: Examples of thermodynamic behavior at various K\, K2 and J3 = 1. 

we found: i) Hint for a first order transition between III and IV in the crossing region of the two 
percolation lines at strong gauge coupling, ii) No transition along the red (vertical) line for small K2 
in agreement with the limit K2 — » 0, iii) Signals for thermodynamic transitions along the remaining 
transitions lines (second order along the horizontal parts of the blue line). 
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4. Summary 

We have observed two transitions associated with the patterns of vortex percolation in the 
Abelian two-Higgs model of high-r c superconductivity. These transition lines are roughly parallel 
to the corresponding hopping parameter axes fQ, almost perfectly parallel at weak gauge coupling 
(j3 = 2). The crossing in some region of finite (JCi, K2) gets new features at strong coupling J3 = 1. 
The Fermi liquid, the spin gap, the superconductor and the strange metallic phases have been 
identified in a wide region of parameters of this model. The percolation transitions are accompanied 
with ordinary phase transitions except at small K"i below the crossing region of the two transition 
lines. First hints for a changing order along the thermodynamic transition lines are found, and the 
joint transition in the crossing region seems to be first order. Further studies are underway to fortify 
our findings. 

More complex Higgs sector will be necessary to build more realistic effective models repro- 
ducing the confinement mechanism of gluodynamics in the sense of Ref. 0. The spin-charge 
separation idea of the high-r c superconductivity should further be extended to gluodynamics to 
reformulate it in terms of condensed matter systems [g|. 
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